home *** CD-ROM | disk | FTP | other *** search
/ BCI NET 2 / BCI NET 2.iso / archives / programming / source / t3dlib_src_r43.lha / writeoff.c < prev    next >
Encoding:
C/C++ Source or Header  |  1995-02-02  |  5.6 KB  |  231 lines
  1. /* writeoff.c - dump the internal database to an OFF file
  2.  *            - written by Glenn M. Lewis - 10/29/91
  3.  */
  4.  
  5. static char rcs_id[] = "$Id: writeoff.c,v 1.11 1995/02/01 22:15:27 glewis Exp glewis $";
  6.  
  7. #include <stdio.h>
  8. #include <ctype.h>
  9. #include "t3dlib.h"
  10. #ifdef __STDC__
  11. #include <stdlib.h>
  12. #include <strings.h>
  13. #include "writeoff_protos.h"
  14. #endif
  15.  
  16. static void process_DESC();
  17. static void process_points_DESC();
  18. static void process_faces_DESC();
  19. static void count_points_faces();
  20.  
  21. static int geom_count;
  22. static char filename[256];
  23. static char *rootname;
  24. static FILE *out, *outcol;
  25. static unsigned long poff;    /* Point index offset */
  26.  
  27. /* Here are a few necessary utilities */
  28.  
  29. static void send_XYZ(f)            /* Print a common string */
  30. XYZ_st *f;
  31. {
  32.     fprintf(out, "%.12g\t%.12g\t%.12g\n", f->x, f->y, f->z);
  33. }
  34.  
  35. static void send_RGB(rgb)            /* Print a common string */
  36. RGB_st *rgb;
  37. {
  38.     fprintf(outcol, "%.12g\t%.12g\t%.12g\n",
  39.         ((double)rgb->r)/255.0,
  40.         ((double)rgb->g)/255.0,
  41.         ((double)rgb->b)/255.0);
  42. }
  43.  
  44. /********************/
  45. /* The MAIN section */
  46. /********************/
  47.  
  48. int geom_flag;
  49.  
  50. int write_OFF(world, name, splitflag, geom_only)
  51. WORLD *world;
  52. char *name;
  53. int splitflag, geom_only;
  54. {
  55.     register OBJECT *o;
  56.     unsigned long total_points, total_faces;
  57.  
  58.     if (!world) return(0);
  59.  
  60.     geom_flag = geom_only;
  61.     geom_count = 0;
  62.     rootname = name;
  63.     if (splitflag) {
  64.         for (o=world->object; o; o=o->next)
  65.             if (!o->extr) process_DESC(o);
  66.         return(1);
  67.     }
  68.  
  69.     strcpy(filename, rootname);
  70.     strcat(filename, ".geom");
  71.     if (!(out=fopen(filename, "w"))) {
  72.         fprintf(stderr, "Can't open '%s' for output.\n", filename);
  73.         return(0);
  74.     }
  75.     if (!geom_flag) {
  76.         strcpy(filename, rootname);
  77.         strcat(filename, ".ipcol");
  78.         if (!(outcol=fopen(filename, "w"))) {
  79.             fprintf(stderr, "Can't open '%s' for output.\n", filename);
  80.             fclose(out);
  81.             return(0);
  82.         }
  83.     }
  84.  
  85.     poff = total_points = total_faces = 0;
  86.  
  87.     for (o=world->object; o; o=o->next)
  88.         count_points_faces(o, &total_points, &total_faces);
  89.     fprintf(out, "%lu\t%lu\t%lu\n", total_points, total_faces, total_points*3L);
  90.     if (!geom_flag) fprintf(outcol, "%lu\n", total_faces);
  91.  
  92.     for (o=world->object; o; o=o->next)
  93.         process_points_DESC(o);
  94.     for (o=world->object; o; o=o->next)
  95.         process_faces_DESC(o);
  96.     
  97.     fclose(out);
  98.     if (!geom_flag) fclose(outcol);
  99.     return(1);
  100. }
  101.  
  102. static void count_points_faces(object, points, faces)
  103. register OBJECT *object;
  104. unsigned long *points, *faces;
  105. {
  106.     register OBJECT *obj;
  107.  
  108.     if (object->extr) return;
  109.     *points += object->desc->pcount;
  110.     *faces  += object->desc->fcount;
  111.  
  112.     for (obj=object->child; obj; obj=obj->next) {
  113.         if (!obj->extr) count_points_faces(obj, points, faces);
  114.     }
  115. }
  116.  
  117. static void process_points_DESC(object)
  118. OBJECT *object;
  119. {
  120.     register int i;
  121.     register OBJECT *obj;
  122.     register DESC *desc = object->desc;
  123.  
  124. /*    if (!desc->pcount || !desc->fcount) return;    */
  125.  
  126.     if (desc)
  127.     for (i=0; i<desc->pcount; i++) {
  128.         send_XYZ(&desc->pnts[i]);
  129.     }
  130.  
  131.     for (obj=object->child; obj; obj=obj->next) {
  132.         if (!obj->extr) process_points_DESC(obj);
  133.     }
  134. }
  135.  
  136. static void process_faces_DESC(object)
  137. OBJECT *object;
  138. {
  139.     register int i;
  140.     register OBJECT *obj;
  141.     register DESC *desc = object->desc;
  142.     register int p1, p2, p3;
  143.  
  144. /*    if (!desc->pcount || !desc->fcount) return;    */
  145.  
  146.     if (desc) {
  147.         for (i=0; i<desc->fcount; i++) {
  148.         /* First check to make sure that this triangle is real */
  149.             p1 = desc->edge[(desc->face[i*3])<<1];
  150.             p2 = desc->edge[((desc->face[i*3])<<1)+1];
  151.             /* if (p1 == p2) continue;    /* How did *this* happen? */
  152.             p3 = desc->edge[(desc->face[i*3+2])<<1];
  153.             if (p1 == p3 || p2 == p3)
  154.                 p3 = desc->edge[((desc->face[i*3+2])<<1)+1];
  155.             /* if (p1 == p3 || p2 == p3) continue; /* How did *this* happen? */
  156.  
  157. #if 0
  158.             fprintf(out, "3\t%lu\t%lu\t%lu\n", poff+p1+1, poff+p2+1, poff+p3+1);
  159. #else
  160. /* Randolf Schultz <rschultz@informatik.uni-rostsock.de> says to do this...
  161.  * Since I do not have a good way to verify this change, I will take his word for it.
  162.  */
  163.             fprintf(out, "3\t%lu\t%lu\t%lu\n", poff+p1, poff+p2, poff+p3);
  164. #endif
  165.             if (!geom_flag) send_RGB((RGB_st*)&desc->clst[i*3]);
  166.         }
  167.         poff += desc->pcount;
  168.     }
  169.  
  170.     for (obj=object->child; obj; obj=obj->next) {
  171.         if (!obj->extr) process_faces_DESC(obj);
  172.     }
  173. }
  174.  
  175. static void process_DESC(object)
  176. OBJECT *object;
  177. {
  178.     register int i;
  179.     register OBJECT *obj;
  180.     register DESC *desc = object->desc;
  181.     register int p1, p2, p3;
  182.  
  183.     /* Process children first */
  184.     for (obj=object->child; obj; obj=obj->next) {
  185.         if (!obj->extr) process_DESC(obj);
  186.     }
  187.  
  188.     if (!desc->pcount || !desc->fcount) return;
  189.  
  190.     geom_count++;
  191.     sprintf(filename, "%s%03d.geom", rootname, geom_count);
  192.     if (!(out = fopen(filename, "w"))) {
  193.         fprintf(stderr, "Can't open '%s' for output.\n", filename);
  194.         return;
  195.     }
  196.     fprintf(out, "; %s\n", filename);
  197.     fprintf(out, "%u\t%u\t%u\n",
  198.         desc->pcount, desc->fcount, desc->pcount*3);
  199.     if (!geom_flag) {
  200.         sprintf(filename, "%s%03d.ipcol", rootname, geom_count);
  201.         if (!(outcol = fopen(filename, "w"))) {
  202.             fprintf(stderr, "Can't open '%s' for output.\n", filename);
  203.             return;
  204.         }
  205.         fprintf(outcol, "; %s\n", filename);
  206.         fprintf(outcol, "%u\n", desc->fcount);
  207.     }
  208.  
  209.     for (i=0; i<desc->pcount; i++) {
  210.         send_XYZ(&desc->pnts[i]);
  211.     }
  212.  
  213.     for (i=0; i<desc->fcount; i++) {
  214.     /* First check to make sure that this triangle is real */
  215.         p1 = desc->edge[(desc->face[i*3])<<1];
  216.         p2 = desc->edge[((desc->face[i*3])<<1)+1];
  217.         /* if (p1 == p2) continue;    /* How did *this* happen? */
  218.         p3 = desc->edge[(desc->face[i*3+2])<<1];
  219.         if (p1 == p3 || p2 == p3)
  220.             p3 = desc->edge[((desc->face[i*3+2])<<1)+1];
  221.         /* if (p1 == p3 || p2 == p3) continue; /* How did *this* happen? */
  222.  
  223.         fprintf(out, "3\t%u\t%u\t%u\n", p1+1, p2+1, p3+1);
  224.         if (!geom_flag) send_RGB((RGB_st*)&desc->clst[i*3]);
  225.     }
  226.  
  227.     fclose(out);
  228.     if (!geom_flag) fclose(outcol);
  229. }
  230.  
  231.